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Quantum algebra symmetry of the ASEP with 
second-class particles 


V. Belitsky and G.M. Schiitz 


Abstract We consider a two-component asymmetric simple exclusion process 
(ASEP) on a finite lattice with reflecting boundary conditions. For this process, 
which is equivalent to the ASEP with second-class particles, we construct the repre¬ 
sentation matrices of the quantum algebra f/^[g[(3)] that commute with the genera¬ 
tor. As a byproduct we prove reversibility and obtain in explicit form the reversible 
measure. A review of the algebraic techniques used in the proofs is given. 


1 Introduction 


The standard asymmetric simple exclusion process (ASEP) ifT^ [T3l |2^ defined 
on a finite lattice with reflecting boundary conditions is a reversible process with 
explicitly known invariant measures. They are not translation invariant, in contrast to 
the non-reversible uniform measures for periodic boundary conditions, and form the 
finite-size analogs of reversible blocking measures (TJ). It has been shown in ifTSl 
that these measures can be constructed by using a non-Abelian symmetry property 
of the generator, viz. its commutativity with the generators of the quantum algebra 

A related process of great interest is the ASEP with second class particles mni. 
For periodic boundary conditions the invariant measures, which are non-reversible, 
can be computed in principle using the matrix product ansatz QEIEI, or alter¬ 
natively using methods from queuing theory a. However, they have a complicated 
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non-uniform structure and no closed-form expression for the stationary probabilities 
as a function of the particle configurations is known. 

The invariant measures for the ASEP with second class particles with reflecting 
boundary have not been studied yet. However, it has been known for a long time 
that the generator of this process has a quantum algebra symmetry, viz. its generator 
commutes with the generators of the quantum algebra t/^[g[(3)] (H]. However, the 
corresponding representation matrices were never computed and the invariant mea¬ 
sures for reflecting boundary conditions, which are expected to be reversible mea¬ 
sures, have remained unknown. In this work we construct the matrix representations 
of the generators of t/^[ 0 [( 3 )] which commute with the generator of the ASEP with 
second-class particles. This approach provides a constructive method to obtain in 
explicit form reversible measures. We also review the algebraic tools required in the 
proofs. 


2 Definitions and notation for the two-component ASEP 


We set the stage and introduce some notation. 


2.1 State space and configurations 

We consider the finite integer lattice A := {1,2,... ,L} of L sites and local occu¬ 
pation variables 77 (k) S S = {0,1,2}. We say that a site k G A is occupied by 
a particle of type A (B) if Tj(k) = 0(2) or that it is empty (represented by the 
symbol 0) if rj{k) = 1. These local occupation variables define the configuration 
TJ = {7?(1),...,T7(L)} G of the particle system. The fact that a site can be occu¬ 
pied by at most one particle of any type is the exclusion principle. 

The following functions of configurations 77 will play a role. Eor 1 < k < L we 
define the cyclic flip operation 

/(J?) = {J7(l),---’?(^-l):’7W + l(mod3),77(k-f 1),...,77(L)}. (1) 

We define 77 *^^ := (/)^^(77) and observe that ( 7 ^)^* = (y^)^. Eor 1 < k <L— 1 We 
also define the local permutation 

a'‘'‘+^iri) = {ri{l),---ri{k-l),ri{k+l),ri{k),ri{k + 2),...,ri{L)}=: 77 “+^ (2) 
We also define local occupation number variables 

■= Sr]{k),0! ^k ^ri(k),li bk'■= ^r](k),2- ( 3 ) 

where 5/t,; is the usual Kronecker symbol with k,l from any set. In particular, we 
define the particle numbers 
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L L L 

N{ri)=Y,ak, M{ri) = Y,bk, V{ri)=Y,Vk- (4) 

k=l k=l k=l 

Notice that N(t 7 ) +M(ri) +M{rj) = L. Occasionally we denote configurations with 
a fixed number N particles of type A and M particles of type B by 

Another useful way to specify a configuration 77 uniquely is by indicating the 
particle positions on the lattice. We write z(rj) = {x,y} with x := {x : ri(x) = 
1}, y := {y : 77(y) = 3}. We call this notation the position representation. Since 
the order of A-particles is conserved we may label them consecutively from left to 
right by 1 to N, and similarly we may label the Z?-particles by 1 to M. By the exclu¬ 
sion principle one has xfly = 0 and by conservation of ordering xi < X 2 < ... < 

1 < yi < y 2 < ■ ■ ■ < In multiple sums over x; and/or y,- such sums will be 

understood as respecting all these exclusion constraints. 0 We note the trivial, but 
frequently used identities 

A(tj) =A(z) = |x|, M( 77 ) =M(z) = |y| (5) 

A^(z) M(z) 

ak= S^„k, I’it = ^ ^,,k- (6) 

i=l i=l 

For a configuration T] = z we also define the number Nkirj) of A-particles to 
the left of a particle at site k and analogously the number A 4 ( 77 ) of B-particles and 
vacancies 14 ( 77 ) to the left of site k 

k-l N(ri]k-l 7-1 M{ri)k-l 

NkiB) := E «<■ = L E M,(77) := E = E E (7) 

!=1 (=1 1=1 i=l 1=1 1=1 

Similarly we define 14 ( 77 ) 77;. 


2.2 The two-component ASEP 

Following 13 the two-component ASEP that we are going to study can be infor¬ 
mally described as follows. Each bond {k,k+i). 1 ^ k < L — 1 carries a clock 
which rings independently of all other clocks after an exponentially distributed 
random time with parameter Xk where Xk = wq if ri{k +1) > ri{k), T^k = if 
77 (k -b 1) < 77 (k) and Xk = 0 ° if 77 (k) = ri(k+ 1) ■ When the clock rings the particle 


* The occupation numbers can be formally regarded as families of mappings >->■ {0,1}, 

b/i :S^ >->■ {0,1} and should thus be understood as functions kj.)?]) of 77 . Since the functional 

argument 77 will always be clear from context (as is the case e.g. in 0), we do not write it explicitly. 
However, we shall usually write explicitly the argument for the particle number functions N{rj), 
M{r]) to contrast them with their numerical values N, M. 

^ We shall use interchangeably the arguments 77 , z, {x, y} for functions of the configurations. When 
the argument is clear from context it may be omitted. 
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occupation variables are interchanged and the clock acquires the parameter corre¬ 
sponding to interchanged variables. Symbolically this process can be presented by 
the table of nearest neighbour particle jumps 

AO -> OA 1 OA AO 1 

OB—5-BO > with rate wg', BO-^OB > with rate . (8) 

AB-s>BA J BA^^AbJ 

We consider reflecting boundary conditions, which means that no jumps from site 1 
to the left and no jumps from site L to the right are allowed. We shall assume par¬ 
tially asymmetric hopping, i.e., 0 < ^ < 0°. By interchanging the role of B-particles 
and vacancies this process turns into the ASEP with second-class particles Ho). 

More precisely, for functions / : C we define this Markov process r\t by 

the generator 

■^/(b) := ^ B')[/(b') -/(b)] (9) 

with the transition rates 

w(b -> b') = L (10) 

k=\ 

defined in terms of the local hopping rates 

w^^^'^{7])=wq{akVk+\ + Vkbk+\ +akbk+\) + wq^'^{'Okak+\ +bkVk+\ +Ma:+i) 

( 11 ) 

for a transition from a configuration 77 to a configuration 77' = 77*^^+^ defined by (|2]). 
The prime at the summation symbol (| 9 ]l indicates the absence of the term 77' = 77 
which is omitted since w(77 —^ 77) is not deflned0 

We fix more notation and summarize some well-known basic facts from the the¬ 
ory of Markov processes. For a probability distribution P{'(]) the expectation of a 
continuous function/(77) is denoted by (/)/>:= Ltj/(b)^(B)- The transposed gen¬ 
erator is deflnedby f{r\) L(j/g§z./(B0-^l7j'(B) where 1 t]'(b) = With 

this definition (| 9 l) yields for a probability distribution P{r\) the master equation 

■^^P{ri) = ^ B)^(b') - w(b ^ b')^(b)]- (12) 

The time-dependent probability distribution P{Tj,t) := Prob [77, = 77 ] follows from 
the semi-group property P{rij) = ^Bo(b) with initial distribution Bo(b) •= 
P{rj, 0 ). An invariant measure is denoted 7 t*{ri) and defined by 

= E^*(b) = T (13) 

n 

A general stationary measure is denoted by n. It satisfies n{r\) = 0 , but no 
assumption on the normalization Tt{r\) is made. 

^ We shall usually omit the set in the summation symbol and simply write . 
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The time-reversed process is defined by 


r>' 


( 14 ) 


with w'^^' irj 77') = w{rj' —?> rj)n{ri')/n{rj). The process is reversible if the rates 
satisfy the detailed balance condition rj') = w{rj' —7 7 ]). We remark that 




( 15 ) 


which is a consequence of (fOl l. 

We define the transition matrix H of the process by the matrix elements 



with w(rj —7 rj') given (fTOl) . In slight abuse of language we shall also call N the 
generator of the process. 


2.3 The quantum algebra Uq[gl{n)] 

The quantum algebra Uq[gl{n)] is the q'-deformed universal enveloping algebra of 
the Lie algebra 0 ((h). This associative algebra over C is generated by i = 
and X^, i= 1— 1 with the relations nmn 


[L;,L ,]=0 
L;X± = 


( 17 ) 

( 18 ) 



( 19 ) 


and, for 1 <ij<n-l. the quadratic and cubic Serre relations 


[X±,X ±]=0 if|/-7|^l, 

(X±) 2 x± -{q + q-^)XfXfXf + X±(X ±)2 = 0 if \i - j\ 


( 20 ) 

1 . ( 21 ) 


Notice the replacement ^ q that we made in the definitions of 0. 



6 


V. Belitsky and G.M. Schiitz 


3 Results 


Before stating the results we introduce for q, q 
number 


q-q-^' 




0 and X €C the symmetric q- 

( 22 ) 


This definition extends straightforwardly to finite-dimensional matrices through the 
Taylor expansion of the exponential. We point out that [—x]^ = —[x]^, [x]^-i = [x]q 
and [x] 1 = X. For integers one has the representation 


A:=0 


Xk—n-\-\ 


(23) 


and the ^'-factorial 



and the q'-multinomial coefficients 


n = 0 
n > I 


(24) 


Cl{N) 


[LW- 

[iV],![L-iV],!’ 


Cl{N,M) 


[LW- 

[N]gm,'-[L-N-M]g\- 


(25) 


The first main result is a symmetry property of the generator. 

Theorem 1. Let H be the transition matrix m of the two-component ASEP defined 
by m with asymmetry parameter q and let , Lj, i = 1,2, j = 1,2,3 be matrices 
with matrix elements 


iYhn’n=^tiiYHk))n'n (27) 

with 

(T2+(^))Vt, = . (29) 

and 77 ^* = ( 77 ) defined by (O. Then: 

(a) The matrices Y^, Lj, i = \ , 2, j = 1,2,3 fonn a representation of the quantum 
algebra t/j0[(3)] (IZZi - (EZD. 

(b) The transition matrix H satisfies [H, Yi^] = [H, Lj] = Ofor i = 1,2, j = 1,2,3. 

The second main result concerns reversibility. 

Theorem 2. The two-component exclusion process rp defined by ^ with asymmetry 
parameter q is reversible with the reversible measure 


( 30 ) 
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Remark 1. (i) In terms of position variables (| 6 ]l we can write 


k{7]) = 


(31) 


(ii) In terms of conjugate occupation numbers = I — aj^^bk = ^ — bk we can use 
the identity ^.[=1 Lf=i (-^jt+i -■*/)= Y!k=\ i2k-L -\)xkX.o write 

n{r]) = . ( 32 ) 


(iii) For finite q one has 
so that n^^(rj) is finite. 


7r(rj)>0 Vrj 


(33) 


4 Tools 

Here we present a review of the tools that are used to prove the theorems. Some of 
these tools are not standard in the context of probability theory. The first subsec¬ 
tion begins with simple facts included for the benefit of readers not familiar with 
the matrix representation of properties of a Markov chain lfT4ll^ . The second sub¬ 
section summarizes more advanced algebraic material from the theory of complete 
integrability of one-dimensional quantum systems. 


4.1 Generator in matrix form 

The defining equation (|9]) is linear and can therefore be written in matrix form using 
the transition matrix (fTbl) 


-^/(J?) = - L (34) 

7)'gS^ 

Notice that the sum includes the term rj' = rj.ln order to write the matrix H explic¬ 
itly one has to choose an concrete basis, i.e., to each configuration rj one assigns 
a canonical basis vector and defines the ordering 1 ( 77 ) of the basis. The set of all 
basis vectors, which we denote by 1 77 ), spans the complex vector space We 
work with a vector space over C rather than over R since in computations one may 
encounter eigenvectors and eigenvalues of H which may be complex since H is in 
general not symmetric. 

Before defining a convenient ordering of the basis we make explicit the relation 
between the matrix representation (fTbl l of the generator and the definition @ of the 
process and rewrite in matrix form some of the Markov properties stated above. 
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4.1.1 Matrix representation of the Markov chain 

Biorthogonal basis, inner product and tensor product: In our convention the basis 
vectors | tj) of dimension d = 3^ aie represented as column vectors. We define also 
the row vector (t? | = | J? ) ^ with the biorthogonality property 

( 7 ? 17 ]') = 5 ^^,. ( 35 ) 

The superscript T on vectors or matrices denotes transposition. 

Consider for arbitrary d two vectors (w | with components w,- € C and | v) with 
components n, S C. We define the inner product by 

d 

{w\v)=Y^WiVi ( 36 ) 

!=1 


without complex conjugation. 

The tensor product |v)(w| = |v)(8)(w| is a d x c/-matrix with matrix elements 
(I v)(w|),,; = ViWj. This notation follows a convenient convention borrowed from 
quantum mechanics. Specifically, we have the representation 

1 = L|J?)(J7| (37) 

of the -dimensional unit matrix. We recall that for two tensor vectors (W | = (wi | C) 
... C) (WL |, IT) = I vi) ® ... (g] IVL) the inner product of factorizes; 

L 

{W\V) = Y\{wk\vk). (38) 

k=\ 

Generator: As a consequence of biorthogonality of the basis one has = 

(7]' \H\ 77 ) and therefore (|9l) takes the form 

■^/(j?) = -(/|7/|j?) (39) 

where the row vector (/1 =J^ri f('n)('n I has components/( tj). A probability mea¬ 
sure Pjrj;) = P{Tj,t) is represented by the column vector 

\P{t)) = Y.P{rid)\ri). (40) 

n 

The semigroup property of the Markov chain is reflected in the time-evolution equa¬ 
tion 

|P(7))=e-"Vo) (41) 

of a probability measure Po{rj) = T’(77o). 

Lowest eigenvalue and eigenvector: Next we define the summation vector 
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(*h=E(j?l (42) 

which is the row vector where all components are equal to 1. Normalization implies 

(s\P{t)) = l Vt>0. (43) 

By taking the time-derivative one has as a consequence 

(sjN = 0 (44) 

which means that the summation vector is a left eigenvector of N with eigenvector 0 . 
This property follows from the fact that a diagonal element of H-q-q is by construction 
the sum of all transition rates that appear with negative sign in the same column rj 
ofH. 

A stationary measure, denoted by | ;r), is a right eigenvector of N with eigenvalue 
0 , i.e., 

//|;r) = 0. (45) 

By the Perron-Frobenius theorem 0 is the eigenvalue of N with the lowest real part. 
The probability vector corresponding to a normalized stationary measure (|43]) is 
denoted by \n*). 

For the stationary distribution we define the diagonal matrices 

Tt* :=Y,n*{ri)\ri){ri\, n ■=J^n{r^)\r^){r^\. (46) 

77 ri 

For ergodic processes with finite state space one has 0 < 7r*(77) < 1 for all 77 . Then 
all powers ( 7 f*)“ exist. In terms of this diagonal matrix we can write the generator 
of the reversed dynamics as 




(47) 


This is the matrix form of ( fTSl l. 

Expectation values: The expectation {f)p of a function /(rj) with respect to a 


probability distribution P( 77 ) is the inner product 


{f)p = {f\P) = {Af\P) 

(48) 

where 



(49) 


n 

is a diagonal matrix with diagonal elements f{r\). Notice that 

/(J?) = (J?I/|J7) = (^I/|J7)- 


(50) 
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4.1.2 The tensor basis 

In order to define a convenient ordering of the basis for H we choose the tensor 
approach advocated in lfT4ll20ll for interacting particle systems. 

Only one site: We begin with the basis for a single site where 77 G §. For a row 
vector of dimension 3 with components w; we write (w| = (wi,W 2 ,W 3 ) and column 
vectors of dimension 3 with components u, we write |v) = (v|^. We define the inner 
product (w|v) := wiUi +w' 2 U 2 + W 3 i) 3 . We choose for a single site the order 11 ( 77 ) = 
1 + 77 and denote the corresponding canonical basis vectors of 

oj, |0) = |2):= l^lj , |B) = |3):= 1^0^ (51) 

With the dual basis (77 1 = [77)^ we have the biorthogonality relation (77I77') = 5 ^ 
The local summation vector is given by (i| := (A| + (0| + (B| = (1,1,1). 

It is useful to introduce the following matrices with the convention | 77 )( 77 '| = 
|77)(8)(77'|; 

fl+:=|A)(0|, fe+:=|B)(0|, c+:= |A)(B|, (52) 

fl-:=|0)(A|, b-:=mB\, c-:= |B)(A|. (53) 

Having in mind the action of these operators to the right on a column vector, we call 
and creation operators, and are called annihilation operators and c* 


are particle exchange operators. 

We also define the projectors 

d:=|A)(A|, 7):=|0)(0|, b:=\B){B\. (54) 

and the three-dimensional unit matrix 

1 = Ll’7)(’7l =“ + ^ + ^- (55) 

n 

They satisfy the following relations; 

a^d = b^d = b^a = c^a = 0, a^d = a^, c^d = c^ (56) 

a^V = b^V = c^V = c^V =0, a+7)=a+, b^v=b^ (57) 

a^b = a^b = b^b = c^b = 0, b^b — b^, c^b = (58) 

da^ = db^ = db^ = dc^ =0, da^ = , dc^ = (59) 

Va^ = db^ = Vc^ = Vc^ = 0, Va^ = vb^ = b^ (60) 

ba^ = ba^ =bb^ =bc^ = 0, bb^=b^, bc^=c^. (61) 


|A)^|1):= 


and 
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With the occupation variables (O for a single site we have the projector properties 


fl|77)=fl|77), 

b\r]) = b\r]), 7)|77) = v|77) 

(62) 

for 77 £ {A,0,B} and a = 13 ( 77 ) = 

Sri^i and so on. Moreover, 


(i|l3^ = {s\V, 

+ 

II 

'ZT' 

+ 

II 

(63) 

= {s\d, 

{s\b = (sj/j, (s|c = {s\d. 

(64) 


L sites: For a configuration 77 £ it is natural and indeed convenient to choose 
the ternary ordering 1 ^( 17 ) = 1 +Lf'=i J7(^)3^^^ of the basis vectors. This corre¬ 
sponds to the tensor basis defined by 


I??) := |t7i)®|i72)®.--® |j?l), (J7|:=(’?iI®(’?2|®---®(i?l| (65) 

spanning the vector space (C^)®^ of dimension 3^. We shall also use the notations 
I z ) and I X, y) instead of 1 77 ). Furthermore, if a configuration 77 has N particles of 
type A and M particles of type B we may denote this fact by writing | r]N,M) for the 
corresponding basis vector. 

The summation vector (i | is given by 

( 66 ) 

This is the row vector of dimension 3^ where all components are equal to 1. The 
summation vector restricted to configurations with a fixed number N of particles of 
type A and M particles of type B is denoted by 

(■sa.m|= 52 (67) 

t]n,m 

The basis vector | rjo Q ) corresponding to the empty lattice is denoted by 10) = 
(■so.o T- 

For matrices M the expression will denote the y'-fold tensor product of M 
withself if 7 > 1. For 7 = 1 we define M®* ;= M and for j — Owe define M®** = 1 
with the c-number 1. For arbitrary 3 x 3-matrices u we define tensor operators 

M^:=l®('^-i)(g)M(g)l®(^-^). ( 68 ) 

Multilinearity of the tensor product yields UkVk+i = < 8 ) [{u (g) 1)(1 ® v)] (g) 

]^ 8 (L-A:-i) _ 1 ^ 0 ( 7 - 1 ) ^ v) and the commutator property UkVi = viUk 

for k^l. 

For u = d or b we note the projector lemma IJl which will be used repeatedly 
below. 

Lemma 1. The tensor occupation operators dk, bk act as projectors 
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N{ri) ^ M(r]) 

akin) =ai,lr]) = £ 7 ]), 77 ) = / 7 ^| 77 ) = ^ 77 ) ( 69 ) 

r=l 1=1 

with the occupation variables a^ and ^ (or particle coordinates Xi and yt re¬ 
spectively) understood as functions of rj or z = rj- 

Remark 2. One obtains the diagonal matrix / ( |49] | of a function f{r\) by the fol¬ 
lowing simple recipe; In f{'(]) one substitutes 77 ,- by the diagonal matrix 77 ,- where 
f\i = at if rii = 0 , fji = Vi if 77 ; = 1 , and fji = hi if 77 ; = 2 . 


4.1.3 Construction of the generator in the tensor basis 


Consider first L = 2 and denote the transition matrix (fTfil l for two sites by h. From 
the definition of the process one computes the off-diagonal part by observing that 
|0A) = |7r'^({A0})) = (a^|A))(8)(a+|0)) = ®a+)(|A) (g) |0)) =aia2\A0) and 

therefore /7 {oa}{ao} = qw{0A |aj"a^|A0) for the transition AO OA, and so on. The 
corresponding diagonal elements h-q-q follow from (fTTI) with the substitution of the 
occupation variables by the respective projectors according to Lemma ([T]i- With the 
tensor basis dfiSl) for two sites one thus obtains the 9 x 9-matrix 


/O 

0 

0 

0 

0 

0 

0 

0 


0 

q 

0 

- q - 

0 

0 

0 

0 

0 

0 

0 

q 

0 

0 

0 

-4' 

0 

0 

0 

-q 

0 

4' 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

q 

0 

-4' 

0 

0 

0 

-q 

0 

0 

0 

4' 

0 

0 

0 

0 

0 

0 

0 

-q 

0 


0 

VO 

0 

0 

0 

0 

0 

0 

0 

0/ 


Next we embed this process on two neighbouring sites in A. By the 

multilinearity of the tensor product the generator becomes /r^tg+i = (g/z 0 

]^8(L-A:-i) Yvjjej-g the hopping matrices 

hk.k+i ■= wq (44+1 - af -f 44+1 “ HK+1 + akh+i - cf c^+i) 

+ wq ^ (4d^+1 - 1 + 44+1 - 1 + bkdk +1 - j X? 1) 

act non-trivially only on the subspace corresponding to sites k^k-\- \ in the ten¬ 
sor space. The off-diagonal elements of hk^k+i are the transition rates h^kk+iq in 
the tensor basis (IfiST l. The diagonal elements of hk^k+i defined in (fTfil l follow from 
probability conservation. The generator for the two-component ASEP on the lattice 
{ 1 ,... ,L} then follows as 

i-i 

H hk,k+i- 
k=l 


( 72 ) 
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As will be seen below, the generator H is closely related to the Hamiltonian 
operator of a quantum spin chain. This is true for other interacting particle systems 
and motivates the terms “Quantum spin techniques” lfT4ll or “quantum Hamiltonian 
formalism” ll^ for the representation of the generator of an interacting particle 
system in the tensor basis. 

Remark 3. Because of multilinearity (l6^ . dMI) are lifted to {s\a^ = (i Itii:, = 
(s|4, {s|c^ = {s\bi, and = {s\dj„ (s|c^ = (sld^.. This yields 

{^\hk,k+i = 0 which implies probability conservation (l44l) . 


4.2 Quantum algebra t/g[gl(3)] and the Perk-Schultz quantum 
chain 

Above we have defined the quantum algebra t/^[g[(3)] ( fT7] i - (l2Tl i. It is convenient 
to work also with the subalgebra t/^[sl(3)]. 


4.2.1 Relation between t/^[g[(n)] and Uq[s{[n)\ 


We introduce generators H,, 1 < i <n and H,- 1 < i < u — 1 through 


H,=H,-H,+i. (73) 

Then the quantum algebra t/^[s[(u)] is the subalgebra generated by and X*, 

i = 1,..., n — 1 with relations ( |2^ . (ISTT i and 


^H,/2^H,/2^^H,/2^H,/2 (75) 

q^'Xfq-^- = q^-JXf (76) 

[X^X7 ] = 57-[H,.],. (77) 

with the unit I and the Cartan matrix 

r 2 if/ = ; 

A,;:= -1 if7 = /±l (78) 

I 0 else. 


of simple Lie algebras of type A„. That Uq[sl{n)] is a subalgebra of Uq[Ql{n)] can be 
seen by noticing that '£!i=i H/ belongs to the center of Uq[Ql{n)] ifTTl . 

We remark that the commutation relations (iTSl l. (l76l l can be substituted by 


[H,-,X±]=±A,7X± 


(79) 

(80) 
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and defining as a formal series through the Taylor expansion of the exponen¬ 

tial. 


4.2.2 Finite-dimensional representations for « = 3 

For « = 3 the quadratic Serre relations (l20l) are void. By introducing 13 

X± := (81) 

the cubic Serre relations reduce to quadratic relations. One has instead of (| 2 TI) 

X± - ^'/ 2 x±x± = 0 (82) 

^i/2x±X± - ^-i/2x±x± = 0. (83) 

and also 

[H,,X±] = ±X±. (84) 

In order to distinguish the matrices corresponding to the three-dimensional fun¬ 
damental representation from the abstract generators we use lower case letters. In 
terms of ( l52b - ( l54l i the three-dimensional fundamental representation of t/g[ 0 (( 3 )] 
is given by: 

X2=b^ (85) 

hi = a, h2 = V, /J3 = h, (86) 

corresponding to 

hi=d—v, h 2 = V—b. (87) 

for the representation of the generators Hi of t/<^[s[(3)]. We also mention the repre¬ 
sentation x^ = 

Next we introduce the coproduct 

4 (Xf) = xf (g) 0 Xf (88) 

4(H;) = H; 0 l-hl 0 Hi. (89) 

By repeatedly applying the coproduct to the fundamental representation, we con¬ 
struct the tensor representations of L(j[s[(3)], denoted by capital letters, 

= Hi='tHi{k) (90) 

k=l k=l 


X;^{k) = q 0X^0 ... 0 , 

Hi{k) = 10... 01 0//,-01... 01. 


with 


(91) 

(92) 
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For the full quantum algebra t/g[g[(3)] we have 

ILL 

Hx = Y.ak=-N, H2='£A=--V, H2='£bk=:M. (93) 

k=l k=l k=l 


Here N and M are the particle number operators satisfying 

^1 = A^l M\riN,M) =M\'nNM)- (94) 

The unit I is represented by the 3^-dimensional unit matrix 1 := 1®^ = N + V +M. 

Notice that Hi (k) = dj^ — Vi; and H 2 {k) = % — bj^. In the L-fold tensor product 
X^{k) {Hi{k)) the term {Hi) is the ^h factor. Therefore (l85l l yields 


Xf{k) = 

X^{k) = 

One has (A'*(A:))^ = 0, X'^{k)Xj (k) = 0 for i ^ j and 

r q^^X^il)X^{k)k<l 
X^{k)X^{l) = \ 0 k = l 

[q^^X^mHk)k>l 

X;^ik)X^{l)=X^il)X^{k)foTi^j. 


(95) 

(96) 


(97) 

(98) 


Thus the spatial order in which particles are created (or annihilated) by applying 
the operators X* (k) gives rise to combinatorial issues when building many-particle 
configurations from the reference state corresponding to the empty lattice. 


4.2.3 The Perk-Schultz quantum spin chain 

We introduce the Perk-Schultz quantum spin chain ifT^ 

L-l 

G=Y.gk.k+i (99) 

A:=l 

where gk.k+i is reminiscent of dlD, but with all non-zero off-diagonal elements 
equal to —1, i.e., 

~8k,k+i = q{akVk+i + Vkbk+i+akbk+i) +9'^' (44+1 +^/t4+i +Mt:+i) 

-^k4+i-l’kl’k+i-^k4+i-4^k+i-l’kH+i-4(^k+i- ( 100 ) 


The gk.k+i satisfy the defining relations of the (3,0)-quotient of the Hecke algebra 
ifTSll which then implies ||6l that 
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[gKk+\.^t] = \gkMuHi]=0. ( 101 ) 

Thus the Perk-Schultz quantum Hamiltonian is symmetric under the action of the 
quantum algebra L(j[s[(3)] and then trivially also under L(j[gl(3)]. 


5 Proofs 


It was pointed out in |[T] that the hopping matrices /td a:+i GD for the ASEP with 
second-class particles satisfy the defining relations of the same (3,0)-quotient of the 
Hecke algebra as the gk.k+\ of the Perk-Schultz quantum chain. This fact implies the 
existence of representation matrices of the generators of L(y[g[(3)] with which the 
hopping matrices /ijt d+i and hence the generator (l72l i commutes. However, in order 
to make this symmetry property useful for probabilistic and physical applications 
one must solve the main problem that was left open in JT], which is to actually 
construct this representation. This is the content of Theorem ([TJ, proved below. It 
turns out that both Theorem ([T]i and Theorem (|2]i are consequences of a proposition 
that we first motivate and then prove. 


5.1 Perk-Schultz chain and ASEP with second-class particles 

The point in case is that G and H differ only by multiplicative factors q and in 
their off-diagonal elements. Therefore the following proposition is a natural working 
hypothesis. 

Proposition 1. Let H and G be the matrices defined in ( 1721 ) and ( 1991 ). There exists a 
diagonal similarity transformation R such that 

G = R^HR (102) 

with an invertible matrix R of dimension 3^. 

Proof. In order to prove this proposition we use the quantum algebra symmetry of 
the Perk-Schultz chain to first construct a good candidate for such a transformation 
and then prove that it satisfies (I1021 i. 

(1) Construction of a candidate R: Fix the numbers N of particles of type A and M 
of particles of type B. By ergodicity the summation vector {sn.m \ dhTl) is the unique 
left eigenvector with eigenvalue 0 of H restricted to configurations with N particles 
of type A and M particles of type B. 

For N = M = 0 one readily verifies that (io,o \H = (■*0,0 |G = 0. The symmetry 
(llOlb then yields (io.o K^i = 0. Therefore, if R exists it has the property 

{SN,M I = YLiN,M){0,0\{Xff{X+fR-^ 


(103) 
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with a normalization factor Yi{N,M). Notice now that (s| I- 

Therefore 


L L-N 

i^\R= L L i^NM\R (104) 

N=0M=Q 

= t L 1l(1V,M)(.o.o |(^rr(^2+r- (105) 

A=0M=0 

Now we make the diagonal ansatz R = R{rj )| J] ) (J] |, normalized such that 

(5o.o|/^|0) = 1 = (^o,oI-R^'|0 ). (106) 

Thus we obtain from (1105b that 

^(ikv.M) =lL(N,M)(.o,o|(^rf (^2+)"I’7 iv.m). (107) 

The normalization factors Yi{N,M) are arbitrary, since they can be absorbed by 
redefining R = RE where £ is a diagonal matrix with matrix elements Yl{N,M) 
in the block N,M. Since both G and H conserve particle number for both species 
one has EHE^^ = El and EGE^^ = G. Therefore G = R '^EIR which implies that 
G = R ^EIR. Hence the Yi{N,M) can indeed be chosen arbitrarily. It turns out to be 
convenient to choose 

Ti(N,M) = ([N],![M],!)-^ (108) 

This reduces the computation of the diagonal elements of R to the computation 

of the matrix elements (0,0 ) from the explicit representation 

( l90b . In order to compute R we first set M = 0 and use 

X 

which is a simple adaptation of an analogous result for the standard single-species 
ASEP taken from dH). The sum over x stands for the sum over all particle positions 
Xi ordered such that x, < xj for i < j, which is the sum over all distinct sets of particle 
positions. Therefore ^^.0 (^,01 = {sn.o \ which allows us to write 

£(x,0) (110) 

and, using ( |69] |, 

(50,01 (Ill) 

Next we apply (^ 2 ^)“ ( |96] | to this vector and observe that for the factors 
that appear in (instead of the that appear in X^) any A-particle is like 

a non-existent site, since H 2 is build from projectors on Z?-particles and vacancies. 
Hence, with regard to the action of (^ 2 ^)^, the existence of A-particles in a config- 
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uration rj^^M = zn,m with N particles of type A behaves like a reduction of system 
size L ^ L = L—N and a coordinate shift 

yi^yi=yi-Ny.{i) ( 112 ) 

for B-particles with Nki^) defined in (|7]i. Therefore the action of {X^)^/[M]q\ on 
(■sq.o I “pvpr yields a ^-factor similar to the one in (II lOl i. but with N replaced by M, 
L replaced by L, Xk replaced by and q replaced by We conclude 

fv-W tv+W 

where the sum is over all distinct coordinate sets and 

R{^nm) = (114) 

Next observe that for any configuration 77 

M{r\) L k-\ L k-\ 

L = E E«'^'t=A^(i?)M(T7)- E 

1=1 A:=U=1 A:=l/=1 

N(j]) 

= iV(77)M(77)- (115) 


For rj = zn,m this yields 

and with Lemma ([T]) 


(■50,0 


(x-f {x+r 

[At],! [M],! 


{snm I?- 


(116) 


(117) 


U=Yi^k-L-\){ak-h) + YY{^i^k+i-b,at+i). (118) 

k=\ k=\l=\ 

Notice that the matrix U does not depend on N and M. Taking the sum over N and 
M then yields from ( 1105b the diagonal candidate matrix 

R = q^^. (119) 

(2) Proof of the transformation property ( 1702b ; We stress that the properties of R 
that we have used in its construction are only necessary conditions for the trans¬ 
formation property (1102b to be valid. In order to prove this property we need two 
more technical ingredients. The first is a transformation lemma, proved in ID 
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Lemma 2. For any finite p f^Owe have 


d‘a^p-^‘ = p^'fljp-^' = af, 

(120) 

<^'bfp-'»‘ = p±^'-h± p“>bfp-^> = bf 

(121) 


(122) 


(123) 


Applying these transformations yields the following auxiliary result. 

Lemma 3. The local creation and annihilation operators transform as follows: 

Ra^R-' = (124) 

Rb^R-^ = ( 125 ) 

RcfR-' = ^±2S=!(^*+i)T2i;Lv+i(4+i)c± (126) 


Proof To prove these identities we use Lemma ^ and commutativity of the pro¬ 
jection operators. This yields 


ELi (aih+i p^ ^k=i TLi ) = 
pLU Tti (±5,,.A+i=F4+i,A)a± = p^ni\ 4±Ekv+i 4)a± 


(127) 


and 

pEti(2<=-'C--i)(4-4)a±^-ELi(2^-^-i){4-4) =p±(2^-2--i)a± (128) 

Joining both yields (I1241 l and a similar computation yields (1125b . Finally, (1126b 
follows from c^ = a^bf, cf = b^a^. □ 

Now we can prove (1102b . We split gk,k+t — stk+i ~ stk+i di^gon^l P^t 

sfk+i off-diagonal part and similarly for hk^k+i- Trivially one has 

^Sk,k+l^ ^ — Sk,k+i — ^t.k+i- 

For the offdiagonal parts, we consider first afaj^^. Eq. (1124b in Lemma (|3ll yields 


RaUf. =fl±^±5(4+i-l)^±2(4-l)flf, 


^k “yt-rr 


k+i- 


(129) 


The general projector property p*™ = 1 -f (p — l)hm. together with (l58b and dMT l 
applied to the subspaces k and A: 4-1 lead to 

(130) 

In the same fashion one proves 

=q^^b^bj^,. (131) 


Finally, by similar arguments 
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Comparing with (fTlT i shows that Rs °;_, *=j and thus ^=H. □ 


5.2 Proof of Theorem ® 

Proof. From (I1241 i. the commutativity of the projectors at different sites, and dSTl i 
lifted to the tensor space, one finds that the local generators Xf[r) transform as 
follows: 


Y+{r) :=RX+{r)R-^ = 

Ff (r) :=RXf{r)R-^ = q-^k=\^'<+^i-r+i^'<a-^ 
Yf{r) :=RX+{r)R-^ = q^'t=\bk-lft^r+i‘>kbf 
Y^{r) :=RX-{r)R-^ = ^-££14+ELr+i 


(133) 

(134) 

(135) 

(136) 


Moreover, since R and Hi are all diagonal one has RHiR ^ ~ Hi. Commutativity of 
the hopping matrices in H with Yf and Lj follows from Proposition ([TJ and 

the symmetry (llOlb of the Perk-Schultz quantum chain. 

To prove the explicit expressions (l26l l. dZTl i for the representations we focus on 
Yf (r). Using the fundamental representation and the factorization property (l3^ of 
the inner product of tensor vectors one finds ( Tj' |a+1 77 ) = ■ The terms in the 

exponential follow trivially from Lemma ([T]i and the definitions Q, (|7]). Following 
similar arguments for the other generators yields the matrix elements of Yf and Lj 
(l26l l. (I 27 ]) as stated in ([T]i. □ 


5.5 Proof of Theorem ^ 

Proof. Since G is symmetric and R is diagonal. Proposition ([T]i implies 


H^ = {RR'^y^HRR^ = R ^HR^. 


(137) 


With ( l47T i we thus have reversibility with a reversible measure k = Rf in the ma¬ 
trix form (l46l) . By the projector lemma ([T]i n yields the reversible measure (l30l l of 
Theorem (|2]i. □ 
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Appendix 

We display some explicit results for unnormalized stationary distributions for small 
lattices L = 2,3,4 and also L arbitrary with small particle numbers A^+M=l,2,3,4. 


At = 0M = 0|00) L = 2 

N=\M = Qq-'^\AQ)+q\OA) 

N^0M= \ q\B0)+q-^\0B) 

N = 2M^Q \AA) 

N=\M= \ q-^\AB)+q\BA) 

N = 0M = 2 \BB) 


N = 0M^0 |000) 

At= 1 M = 0^-2|^00) + |OAO)+^2 |ooa) 

N = 0M=1 q-^\00B) + \0B0) +q^\B00) 

N = 2M = 0 q-^\AA0) + \A0A)+q^\0AA) 

N=lM=l q-^\A0B)+q-^{\AB0) + \0AB))+q{\BA0) + \0BA))+q^\B0A) 
N = 0 M = 2 q-^\0BB) + \B0B) +q'^\BB0) 

N = 3M = Q \AAA) 

N = 2M= \ q-^\AAB) + \ABA)+q^BAA) 

N= \ M = 2 q-^\ABB) + \BAB) + q^\BBA) 

N = 0M = 3 \BBB) 


N = 0M = 0\0000) L = 4 

At= 1 M = 0^-^|A000)+^-‘|0A00)+^|00A0)+^^|000A) 

N = 0M=1 q-^\000B)+q-^\00B0)+q\0B00)+q^\B000) 

At = 2 M = 0 AAOO) + ^-2|^0A0) + IAOOA) + IOAAO) + oaoA) + ooAA) 

At= 1 M = 1 ^-5j^OOB)+^-3(|AOBO) + |OAOB))+^-'(|ABOO) + |OABO) + |OOAB)) 
+^(|BA00) + |0BA0) + |00BA))+^3(|50A0) + |0B0A))+^5|B00A) 

N = 0 M = 2 q-^^lOOBB) +q-^\0B0B) + \B00B) + \0BB0) +q^\B0B0) +q‘^\BB00) 

N = 3 M = 0 q-^\AAA0)+q-'\AA0A)+q\A0AA)+q^\0AAA) 

N = 2M= \ q-^\AAOB)+q-^{\AABQ) + |A0AB)) +^-'(|0AAB) + |ABA0) + |A0BA)) 
+^(| ABOA) + IOABA) + | BAAO)) + BAOA) + | OBAA)) + ^^|B0AA) 

At = 1 M = 2 ^-5|A0BB)+^-3('|^50B) + |0ABB))+^-'(|BA0B) + |0BAB) + |ABB0)) 
+^(|0BBA) + |BAB0) + |B0AB))+^^(|BBA0) + |B0BA))+^5|BB0A) 

At = 0M = 3 ^-3|0BBB)+^-'|B0BB)+^|BB0B)+^^|BBB0) 

At = 4M = 0 |AAAA) 

N = 3M= \ q-^\AAAB)+q-^\AABA)+q\ABAA)+q^\BAAA) 

N = 2 M = 2 q-^\AABB) + q-^\ABAB) + \ABBA) + \BAAB) + q^\BABA) + q'^\BBAA) 
N = \ M = 3 q-^\ABBB) + q-^\BABB) + q\BBAB) + q^\BBBA) 

N = 0M = A \BBBB) 
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N + M=\\ 

;r*({x},0) oc ^2-1 

;r*(0,{y})oc^-2,+i 

N + M = 2\ 


7r*({xi,X2},0) 2 

r «2xi-2yi-l 

^*({^i},{yi}) 


2 xi- 2 yi + l 




yi < 2 ci 

yi >xi 


N + M = 3: 


7t*{{xuX2,X3},9) ^ 


'^*{{x\,x2},{y\)) 


‘^*{^,{y\^y2,y2}) 


^2xi+2x2-2_V|-3 

^ 2 xi+ 2 x 2 - 2 v |-1 

2xi+2x2-2vi + 1 

H 

^ 2 x 1 - 2 yi- 2 ^ 2-1 
^ 2 x 1 - 2 yi- 2 ^ 2+1 
^2xi-2yi-2y2+3 

-2yi-2y2-2y3+3 


yi <xi,X 2 
xi < yi < X 2 

2 ti,X 2 <yi 

yi,y 2 <xi 
yi<xi< y 2 
XI <yi,y 2 


N + M = 4: 


a:*({xi,X 2 ,X 3 ,X 4 }, 0 ) oc 


q 


'2,x\ -1-2x2 “I“2 a'3 “1~2x4—4 


7r*({xi,x2,X3},{yi}) oc 


;r*({xi,X 2 },{yi,y 2 }) 




{ ^2x1 +2 x2+2x 3 -2yi -5 
^2x1 +2 x2+2x 3 -2yi -3 

q2xi +2x2+2x 3 -2>'i -1 

^2x^ +2 x2 +2 x3 —2>4 +1 
' ^2xi+2x2-2yi-2y2-4 
^2x1+2x2-2yi-2^2-2 
^2xi+2x2-2yi-2y2 
^ ^2xi+2x2-2yi-2y2 
^2xi+2x2-2yi-2y2+2 
^2xi+2x2-2yi-2y2+4 

q2xi -2vi -1 

^2xi-2yi-2y2-2y3 + l 

^2xi-2yi-2y2-2y3+3 
^2x1 -2yi-2y2-2y3+5 

2yi -2y2 -2y3 -2y4+4 


yi <xi,X2,X3 
XI <yi <X2,X3 
xi,X 2 <yi <X3 
X1,X2,X3 <yi 

yi,y2 <2ci,2C2 
yi < xi < y 2 < 2 C 2 
yi <xuX 2 <y 2 
xi <yi,y2 <X2 
Xi <yi <X 2 < y2 
XI,X2 <yi,y2 

yi,y2,y3 <xi 

yi,y2 <xi <X3 

yi <xi <y2,y3 
XI <yi,y2,y3 


T^*{^,{yi,y2,y3,yA}) q 
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